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ABSTRACT: 

As S0 i.s are the most lively candidate for the so.ai.ed missing carbon sink, 
we explore the possible impact of C02 fertilization on the global humus inventory. 

For 'inv °iven greening-induced enhancement of plant growth, 

^i^wi.ldependonthespecmumoftumoverdmesMthrespectm 


oxidation. Here, we 

measurements. 


develop estimates of carbon turnover rates 


from soil radiocarbon 


INTRODUCTION: 

The so-called missing carbon sink is the difference between the amount of 

carbon released by fossil fuel burning and deforestation, and the. amount 

carbon which has appeared in the atmosphere and ocean reservotrs tgtn 

„ When combined with the atmospheric increase, most estimates o - 

uptake (Keeling, 1973; Oeschger et al„ 1975; Broecker et al„ 1979; Bacastow and 

,981- Pena 1986; Maier-Reimer and Hasselmann. 1987; Keeltng et ., 
Bjorkstrom. 1981, Pen . shon of matching fossil fuel C0 2 

1989a b and Toggweiler et al., ) 

• ie interval Thus, the missing sink is generally thought to 

released during the same time inters . 
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r, a hle 11 If Houghton's (1990) estimates of 
be st leas, as large as deforestation (table 1). If ^ ^ ^ ^ 

-in* used (figure 1 and table 1). ewer P 

biospheric releas 1992) are unaccounted for. This 

5 rrr/vear Surmiento and Sundquist, 

37 G.C (or 1.8 GtC/year, fenUiz ation and N depostuon, 

bookkeeping does not include processes such as C02 

Which could stimulate terrestrial carbon storage. of „ e terIes mal 

One promising prospect for the mtsstng si phere contains 

. ,, h m, fertilization of plant growth. Todays a 

biosphere through CO. in 1800 (Neftal e. al„ 1985 and Bacastow 

25 % more C02 than did the atmosp ere ^ increase 

and Keeling. 1981). 

. divided by percent increase in atmospheric v. 
in plant growth r s 0.35 (Esser, 1987. Kohlmaier e. ah. 1989 and 

short term chamber expel intents 1 lhe natural environment, 

1 locm If these experiments are apput 

then between 180^ hotvterm chamber experiments to estimate real-world 
BUt ' ' 

^b, nitrogen availability, htdeed, the long-term experiment a, 
limitations impose „„ signifto „t increase in pian. biomass for 

msuits of Norby « al. - . ^ envir o„ments; however, these 

yellow-poplar trees grown lor . dioxide ^ ^ rf 

f increased root turnover and increaseu ^ 
authors do repor ... least for some regions. N is no longer 

soil, deposition of NH 3 andN03 (Kauppi e, al., 1992). 

limiting cst jmate that 0.05 Gt of oxidized N and 0.054 Gt of reduced N are 

' , j ip the oceans annually. Even if a relatively low 

deposited on lan , ed utilization of this nitrogen by the global 

carbon/nitrogen ratio of lOts adopted, utthza, ton 

ecosystem would sequester 1G.C per year. 


MODELING APPROACH: 


Part0 „ « al. (1987 and 1988) have developed the Century model and 
Jenkinson (1990) has developed the Rolhemsted model for studying soil 
matter inventories. Die structure of the two models ts similar, or examp e. 
estimate the flux of p.ant residue into three sod compartments: acttve. slow an 

passive.^ C(;ntur y mo del calculates the climate change response of soil organic 
matter invents for a S rassland ecosystem, .fs no, clear tha, these results can he 
applied to other ecosystems. The parameters showing the greatest sensitivity are the^ 

indirectly for long-term incubation experiments. 

The Roihamsied model’s structure is similar to the Centun mo 

Jenkinson estimates his turnover times by using data from carbon inventory changes 

in lone-term agricultural expenments : 

Here we present a strategy for es.inxt.ing turnover times more directly us, 

radiocarbon measurements. To explore the dynamics of COa fertilization’s influence 

on son carbon siorage. we have deveioped a simple mode, We assume that - 

carbon inventory for the world's soils was a. steady state in the year 1850 Ox, 

input flux via litter and roots etjuaied the output flux via bacteria, decomposition,. 

perturb the system by increasing the level of atmospheric CO^ 

time history from 1850 to the present) and assume a C02 fertilization factor 

, T> .: s in lurn causes the inventory of soil carbon 

for the carbon flux into the soil box. This, in turn. 

,o rise However, as the carbon inventory increases, so does the rate o sol 
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C0 2 in each year's time step). We further assume ^ the change tn glo 
temperature over this interval has been too small to increase rates. So, 1 carbon 
accumuia.es because the rise in the rate of decomposition lags the me ^ 
input This lag depends on the spectrum of residence times of car n atom 

Clearly, no single residence rime can characterize the complex of organic 
compounds in soils. Nor can the most appropriate residence time for wry ^ ^ 
applied to others. Even the determination of the portion of the soil carbon ^allmS mo 
broad residence time bands currently lies beyond our grasp. As outltned . 
simplify the problem by dividing soil organics into that porrion which has a mean 

j „ „„ m — •• ' - “I — *«• 1 - “ ■ “ 

^‘^urbntions. like COJ fertilization, have the gma.es. impact on carbon pools 

having short residence times, and the leas, influence on carbon pools havmg long 

residence rimes. For example, .he carbon poo, having n residence „me of a ew^ears 

increases its carbon invenrory by 9, (figure 2, in conmsr. poo.s wrrb res, dence 

times greater than ,000 years experience almost no change. A carbon pod havmg a 

25-year residence rime will show an increase of 5.3%. 

Tracer experiments provide a way » esdmare soU-carbon residence nmes. 

Two giobal carbon tracer experiments include diluting and enhancing the C-14 C °^ m 
of the atmosphere. We define R as the measured C- 14/C ratio of the 
by lhe c ,4/c ratio of 1850 wood. R(armospbere) debased by 15% berareen 85 
^ ,950 as rhe result of burning fossii fuel, this diiurion's impae, would be drificui, 
to detect However, nuclear bombs increased R(soii) from 1950 to 1963. to almost 
, 9 In ,96a. nuclear rests were banned and since rhen Riarmospbere, htrsdecreased 
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ommdally This near-doubling of R(atmosphere) should be reflected in the fast 


pools have ever 


smaller maxima occurring at ever later times. 


EXPERIMENTAL RESULTS: 


We have obtained radiocarbon data for native sites from four locations. . 

Italv . Duke Experimental Pores, neat Outham, North Catolina; the Konaa *** - ^ 

' . . Kansas' and two sites near Saskatoon, Saskatchewan, Canada (see ta e- 
Manhattan. •• picking of visible plant 

rt c ch oenau 1988). Our sample treatment included hand pu: 

I eln removed with dilute acid. The remain.ng orotic carbon has been ox, diced 

rlnd^eandsubsepuently reduced, o form graph, ,ar g e G ,r acceleramr 

a-—— 

t'.-’ ““ b “ 

rensonablv comprehensive, data are scant. 

When the native soil radioed values are plotted against time (figure :,), a 
, UniBd number of pre-nuclear soils show values below 

“"d bet0n ' iS JPP '' rcn 0 () ( ,, lble ,, Tlte values rise during die 1960's and flanen 
atmospheric averaging 0.90 (table 

thereafter at a value of 1.15. 


^ t *„ sv cv ::*J k H ; 


interpreting radiocarbon measurements using model 
RESULTS: 

» — - — “ ' 1 „ 

. present-dav values (figure 5). For example, the average pre 

le (table 2, is 0.90 (850 years,. Ourmode! predicts thatearbon having atumover 

thnTof 850 years will show little bomb C-M build-up, but R(soil)observations 

I! , , Hence, sol, carbon poois must have a continuum of exchange rates 

l uivelv fast We have modeled the soil carbon system 
ranging from very slow to relatively fast. We 

one slow and one last. 

. f . , s s0 as to fit the mean of the observations. our first attempt to do this, 
,hese radio Qf 3 J00 yC ars to the slow fraction and 25 

W e toe :t:' ^ f rvemahe,hefa,frac,ion 75 ,of,he,o.a,and,heslow 

POSt ' b CetreTirbre Jown as follows: We assume dia.deep in the soil zone 

uenre .he R measured deep in the soils provides 

,he fast-cvclina component is gone. Hence, the Kmeu 

the fast „ nme The available deep soil 

nf , he slow-cvcling component’s turnover nme. The aval 

T’T on values listed in table 3. average 0.63, which corresponds to a 3.700-year 
residence time. For the fast-cycling pool, we find by iteration that a 25-year residence 
time pool fives the best fit. Specifically, we show how bomb radiocarbon won, 

time pool s , 00 . 25- and 10 -year residence nmes 

build up in soil carbon pools having 3’^^’ 75 . 31 , p^.he fast pool to replicate the 0.90 

(figure 6 ). We add 25% of the slow pool observations over our 

pre-bomb R(soil, value. We supenmpose -he expenm 
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■ find that a 25- year residence time for the ^ P 00 * provides the 
model prediction to find that a 4a y 

best fit to the data (figure 7). dividi „g , he fas, turnover 

A better fit could prefab y => ^ ^ ^ , im , Wc show one such 
reservoir into several sub-parts, eac « ^ ^ ^ ^ and the other 

example in figure 8. A combination ^ significant improvement. Our 

with a 100 -year turnover ume oes no is „ ol warranted until a far 

conclusion is thatfurther subdivision of the fas, reserv 

larger data set IS available. . advance and provide a template for future 

While our approach is crude, r, oes^ ^ ^ ^ soik and bo mb-era 
research. If a larger set of measurements on fw speciflc 

soils were available, we could make esrimt»te8^| collected in the early I970‘s, 

USING CARBON RESIDENCE TIMES TO ASSESS THE IMPACT OF C02 

FERTILIZATION: . md residencc ti me. we can estimate the 

Based on the rcd because of COa fertilization, 

amount of additional carbon that carbon, we need to know the fraction 

TO calculate the global " ^ Usjng soil rad iocarbo„ and carbon 

of giobal sod contain, ng rap, ^ ^ ^ ^ cycling carbon, 

content profiles, we esumate th ^ „ ,* We 

Schlesinger 0991) estimates that ^ Bms wi , h small amount of 

subtract 500 G,C from ^ and' swamp ecosystems (table 4). 
fast -cycling ettrbon, ^ ^ cMbon pool » be 450 OtC (1000 G,C 

rimes O^S.^ahlc 6). Using a 75 -year residence time and applying a 0.35 greening 




w Vi*»*>*< 


• i «fO for 7 GtC from 1958 to 
factor, this translates into a 24 GtC seqestrauon since C 

1978 or 0.6 GtC/year during the 1980’s). 

This approach ignores litter and fine root storage. In our study, ese 

components have been physically temoved front the soli. 1* »~P«* ° 

55 GtC (Schlesinger, 1991) and has an average residence time of three years (Esser, 

i987). A three-year residence ante translates into a 7% carbon 
(able 5). If we assume about the same input from fine root turnover, 
fine roots and liner have been sequestered because of C02 fertilization Sln ° C J^ ^ 
We have also failed to consider the increase in vegetation btomass. 
a , 5 .year residence time for short-lived vegetation having a 90 GtC inventory an a 
63-year residence time for a long-lived 500 GtC pool (Wameck, 1988), WC 
that the short-lived pool sequesters 0.2 and the long-lived pool 

annually (figure 9). ,f we combine the amount of carbon sequestere 

fertilization in soil, liner, fine mots and vegetation, we calculate that about 5 GtC 
have been stored from 11150 to 1985, 17 GtC from 1958 to !978, and 1.5 
annually for the 80’s (table 5). Themfote. about 80% of the annua, mtsstng 
1.8 GtC/year (IPCC, 1990) might be explained by C0 2 temhzauon. 


CONCLUSION: 

We have developed a mode, that estimates soil carbon residence times from 

soil radiocarbon measurements. We use these residence times to ^ 

soil carbon's response to C02 fertilization. We hope dtis appmach wt.l . 

strategy for future research. 
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Table 1: Global carbon budget. 

sources. 

gtC Eeference 

fossil fuel 76 Marland & Rotty, 1984 


amBllidsLiaffi 
Gic EsIeisdss 

5.4±0.5 


IPCC, 1990 & 
1992 


deforestation 

32 

Houghton et al, 1990 

l.6±l.O 

7.0 

total additions 

108 



sinks 


Bolin, 1986 

3.2±0.1 

atmosphere 

43 

increase 
ocean uptake 

28 

Broeckeretal 1979 

2.0±0.8 

total storage, 
missing sink 

71 _ 
37 

(by difference) 

1.8±0.4 


Table t Soil Radiocarbon data for non-cultivated soils. 

R(soil) location 


sampling 

jfolf- 


1927 

1959 

1956 

1959 

1960 
1962 * 


1963 

1964 

1965 

1965 

1966 

1967 

1968 
1969 * 


1970 

1970 

1970 

1970 

1971 
1971 
1971 
1971 
1973 
1977 
1977 
1984* 


1984* 

1985 

1991 

1991 


0.82 

0.90 

0.88 

0.93 

0.95 

0.89 

JLS6- 


Leningrud, USSR 
Amazon basin 
New Zealand 
Amador County, C A 
Judgeford, New Zealand 
Saskatchewan, Canada 


vegetation/ 
SOiLUOK — 


podzol 
ultisol 
silt loam 
temperate forest 
silt loam 

cherozemic 


1.01 

1.02 

1.00 

1.05 

1.07 


Judgeford, New Zealand 

it 


silt loam 


precip. 

estimated 

fmm/yr) £m-°C- 

325 l n 

1500 30 

1250 15 

250 

1250 15 


325 


1250 


15 


1.08 

if 

R 

»• 

*7 

1.11 

1.01 

n 

Saskatchewan, Canada 

cherozemic 

H 

325 

tt 

7 

n 

m 

0.99 


It 

ft 

»t 

1.12 


ft 

If 

H 

1.07 

11 


It 

If 

1.04 




l« 

1.06 


M 

It 


1.12 


M 


n 

1.03 


It 



1.03 


11 

If 

M 

1.12 


ft 



1.13 

11 

ft 

" 

H 

1.09 

It 

II 


8 

If 

1.19 

1.13 

1.12 

1.12 

t.ll 

1.11 

1.12 

1.09 

1 -IT 

it 

Wohldorf, FRG 

hapludalf 

H 

650 

n 

Ohlcndorf, FRG 
Timmcndorf, FRG 
Akka. Israel 

hapludalf 

polloxcvert 

620 

IT 

•f 

18 

tt 

Qcdcma, Israel 
Patanchcru. India 
Lodi. Italy ’ 
Saskatchewan, Canada 

rhodustalf 

grassland 

grassland 

grassland 

deciduous 

1300 

600 

325 

732 

27 

14 

7 

13 

u5 

1.12 

Konza Praric, Kansas 
Durham, NC 

1100 

15 


f 

f 

a 

n 

l 

e 

_S 


1 

1 

d 

d 

i 

i 

i 


i 

j 

j 

j 

1 


1 

g 

g 

b 

b 

b 

b 

b 

b 

m 

m 

m 

m 


f| , . bv an a „ , he act „al dates the soil was collected. Those sampling dn.es followed by n 

Dates followed by an ^ ^ years from the publication date. We define pre-bomb as being earli 

have been estimated. d(M artel & Paul. 1974b); 

than 1963 . Refer el a!., 1976>; j(Goh e. a 

„ u it Af ,1 1967V fCTrumbore et al., 1990), g(Ounen, tvo ,, 
e(Campbell * , 1997 ') In Kansas, we sampled the 

19 77 ); KO Brien & Stout, 1978 ); mOhis study) and n(Trumbore. 199 .). 

Konza Prairie, an NSF LTER site. 



Table 3: Deep soil radiocarbon values. 

soil type 

deElhlsul 

jUsoiD 

Mollisol 

65-105 

0.49 

Mollisol* 

40-50 

0.76 

Vertisol 

60-140 

0.60 

Udic 

60-140 

0.60 

Mollisol 

85-110 

0.62 

Mollisol 

80-100 

0.71 

♦under house. 



China 


Beckcr-Hiedmannetal, 1987 


New Zealand O’Brien, 1986 

Israel Scharpenseel and Becker- 

fcrael Heidmann, 1989 


Patancheru, Bccker-Htidmann, 1989 
India 


Hildesheim, 

FRG 


Tsutsuki et al., 1987 

n 


Tabl e 4: Soil organic carton inventories, G,C (S< =“" ^oi, carton 

accumulating ecosystems, whose carton turnover ,s very slow f 
inventory to gain a more accurate estimate of the global fast carton pool. 


Reservoir 

total soil organic carbon 
inactiveanls^ alpine , 5Wa mp) 

active soils 

fast pool 
slow pool 


Tarhon inven tory (GtQ 

1500 

500 

1000 

450 


550 


Table 5: Carbon inventory increases. 

residence 
time 

ffitO— _*J£2TSL 


carbon pool inventory 
size 


fast soil 

litter, fine 
roots 

long-lived 

vegetation 

short-lived 

vegetations 


450 

110 

500 

90 


25 

3 

63 

1.5 


1850-1985 


5.3 

7 


3.5 

8.5 


total 

missing sink 

* • e-ini- Hprreases to 23 GtC from 1958 to 
*The missing sink decreases - 

estimates of biosphenc release are used. 


*8^ .958-1978 — 


24 

7 

18 

8 


57 

105 


7 

3 

5 

2 


17 

23*-37 


IGiCZvrl 

0.6 

0.3 

0.4 

0.2 


1.5 

1.8 


1978 if Siegenthaler and Oeschger's (1987) 



Figure 1, Tide: Missing cartonsink, 1958-1978 

, j-r nresta iion) and sinks (ocean uptake, atmosphenc 
Sources (fossil fuel an ^ kft ^ 

- r -utmrwoheric caibon dioxide from l^ s 10 17 
increase and "missing) for atmospn fossil fuel and 

deforestation exceed ocean uptake and the aflnosi) “ WOTtbeaI bonin 

,, 5 onm in 1958 to only 335 ppm in l?/o- w 

atmosphere increas ^ me th e C 02 content in 1978 would have been 352. 

the missing sink .0 have remained anbome. 


CO 2 (ppm) 





Figure 2: 


Title: Soil carbon residence time vs percen 


t carbon reservoir increase, 1850 to 1985 


This plo, she*, how residence dm. influence, d* aA “' 

inventoiy fora 0.35 C02 fertilization factor. The shortest cycling time caibon pools have 


the greatest increase. 



C/C/ 14 C^1850 W ood 




Figure 3: Title: Mootbo" specmun pl°t 


, ,aca ( 0 1995. The ratio ranges front 0.8 to 

C44 redo is ''0"^“ rgreu.es, increase, .Hch stems fan nuclear 
^,os, 1.9. ^e^ ^^b^e^p^^^ghUtesoil^ltavuig 

testing in Ae atmosphere. soil caibon residence time plotted, 3 years, shows 

different ^^k****^^ However, even in this shot* residence time, the atmospheric 
the greatest pettu ^ nc of \ 000 yean shows the least 

signal is significantly attenuated- Tiie longest residence tun 









R(soil) 



Figure 5: Tide: R(soil)vstiine 


original 

Of poor 


FAGL' i $ 

QUALITY 


Measured R(soiI) caitxxi values for non-cultivated soils ait ploned against time 

f m „50,o 1991. The values tend 1° increase during the , 9 ® Btrieve, off. 

^ . ocq vear residence time are shown by 

Model results for a theoretical aAon P°° leaving an 850-year resi 

the line. If soil caibon consisted of components having this single residence ume, one 

. ,. lf ; ,. ns to ” n’DunH this 



Rsoil 


. v V • 4 ii.vT*vil 



Figure 6: Tide: toA taving residence,^ of 3700, , 00 , 25 and ,0 years vs rune 

Soil radiocarbon values can be modeled as a system consisting of slow and fast 

— — 


R(soil) 



h . nc 75% of to-cycling components wiO, 25% of the slow^cl^S 

. curves diverge when bomb radiocaibon is released to the armosphere. 

However, these curves aiveig Wehave 

fnmon-cuitivated soils is also plotted (open 

R(soil) caibon . . 5 . 25- and 40-year residence times 

• rw&A our soil carbon modeling results for 15 , •« an 
superimposed ^ ^ ^ ^ tns , bracket :ne experim tntal data . 


converge because 1 



t ( / 


:^\}hr++ 


»A W-V 



ure 8: Title: 
dictions 


5 Comparison of single fast pool wilh a twoa>rnponent fas' pool 


. „ „f *= model calculation forour aamlari tw^mpoacnt 

nsisting of 25% 3700-year, 45% .O-yearand^ l«H- 


percent carbon reservoir increase 



Figure 9: Title: 


ft*, carbon resetvoir increase forvariouscaAonpools. 


We plot the pen^t reservoir increase from i^OO t° 1990 for atmospheric c °2 

observations, short-lived vegetation, fast soil and long-lived vegetation. 
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AS first recognized by Heinrich (1988), a series of deposits rich in ice-rafted 
debris and unusually poor in foraminifera are present in northern Atlantic 
sediments of the last glacial period. Subsequent studies reveal that these 
deposits form a belt extending across the entire Atlantic Within four of die 
deposits are layers with an unusual set of properties. Each was deposited 
rapidly, each contains abundant detrital carbonate, and each has day-sized 
material with unusually high K/Ar ages. The evidence leads us to propose 
that each layer records massive discharge of icebergs originating in Canada 
and possibly northwestern Greenland. The deluges appear to be the product 
Of sudden advances of ice streams onto continental shelves combined with 
cooling of North Atlantic surface water. The catastrophic events may have 
been a consequence of dramatic cooling events or stochastic surging of the 
Laurentide-northem Canadian ice sheets. 


. , ... , ludv of the sediment composing the 
To date the most detailed study 

^ ”„ o r, .. 

— r * — - - 7'“ ;r;r; z~ •* 

c loarlv define the six Heinrich deposits of the last gi Y 

-21,000 and -28,000 years), and by extrapolah g ^ and 

we have estimated the ages for the remaining 3 laye , - 

-60,000 years; Fig. D- AMS radiocarbon ages nearly ^ B J 

609 have been obtained for layer 
7 c 55 f rom the Labrador Sea (Fig. 1; Table 1). 

At DSDP site 609 four of the Heinrich deposits (1, 2, 4 and 5) con am 

th1 v ravs than the adjacent sediment, 
prominent layer that is less transparent to X-ray 

, V « „ lithic fraction in these four layers reveaiea 
Grain counts of the > 6 P . jgy (Fie. 

^singly large amount of detrital — ^ _ 

»• ^ detrita ‘ Calb0nate 8tamS “ few contain oolites and fragments of 
grained recrystallized textures. A ew „ that 

echinodems. Petrographic examination the P 

from 40% to 60% of the identifiable matenal is also detrital 

f 40 • t the distinctive appearance of the layers in 

abundance probably accounting for the distinctive pp 

X radiographs. In contrast, the mineralogicai composition of the other two 

8 P • m that of the ambient sediment, consisting of 

T-Tpinrich deposits is similar to that ot tne * 

j feldspar variable amounts of black volcanic g ass, 

60% to 90% quartz and feldspar, va 

a locc than 8% detntal carbonate (rig. *■)• 

-° r —nrr;: "lid ^ ~ - — 

. . , / 1988 ) ij n e counts of the > 63 micron fracti 

seamount cores studied by Heinrich (1988). Line 
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reveal the same unusual abundance of dehital carbonate in layers 1.2,4 and 5 ^ at 

DSDP site 609 (Fig- 2). These layers also have unusuaily >° W p«os^ ^ ^ 

about one half that of typical deep sea muds. Perhaps e m0 ^ ^ 

these cores is that m all four layers witn + g 

average K/ Ar age of 897 ± 15 Ma, contrasting sharply with average ages of - 

\/fa fnr the < 2u ambient sediment (Fig. 2). 

r ,;j ru , - ».— 'V- - * — - “• -rz, 

between about 45°N and 52»N (Fig. 1), essentially the same area o >g 
accumulation previously identified by Ruddiman (1977). W.thm th.s e 
thicknesses of the carbonate-rich layers increase markedly to the w s ^ 
cores from the Labrador Sea. Andrews and Tedesco (in rev.ew) and T n 

Andrews (1992), documented even thicker equivalents of the carbonate ®D 

• , , ,, *»t the source of the detrital carbonate must 

layers 1 and 2, and they concluded that the 
y uoHrnrk in Hudson Strait (Fig- 2). bimu 

have been nearby limestone and dolormte bedrock m n 

carbonate rocks are widespread in eastern Canada and northwest Green an ■ - 

nd they overlie a crystalline basement with ages old enough to be the source o 
Z K Ar ages hr the carbonate io in - open ocean. Thus, the Labrador , 
It have been one of the principal routes through which icebergs carrymg thrs 

distinctive sediment entered the North Atlantic. 

Two other features of the Heinrich deposits are worthy of menhom 
Pics., the layers with abundant carbonate 1RD must have bcen depos.te 
rapidly relative to the ambient sediment. X-radiographs of cores from DSDP 
site 609 and the Dreizack area indicate that each layer has a sharp base m 
distinct contra, to the gradational bases of all other «D layers observed to the 
depth of the 5/6 boundary. In addition, the flux of l.th.c grams 
carbonate 1RD layer in H2 (the only detrital carbonate layer in whrch we 


3 



make a flux estimate so far), was much higher than that of the sediment on 
make a nux « evidence of dissolution of 

either side (Table 2). Although we have fo 11- 

forams in H3 and H6, forams in the iayers with carbon te 
preserved, and dissoiution does not affect our estimates of sheil ram 

faUnal S «In P d. all six of the Heinrich deposits formed during extreme climatic 

JZ ; OSOP Site - and in the Oreizack area (Heinrich, ,SS the 

low foraminiferal layers are dominated by especiaily large abundances of the 

, « \ (p; e 2‘ Heinrich, 1988), indicating 

P0 ' ar U Vema core ^ - the 

ncntheastern Atiantic (Fig. a, stages 2 through 4 contain 

W ith nearly 100% N. and nearly 0% coccoht s, m^a ^ 

extreme southward penetration of polar water (McIntyre a ^ 

, d that three of these intervals contain carbonate an 

found that three other three (McInty re et 

correlative with H2, 4 and 5. Approximate ag 

, c nf HI 3 and 6. Our estimates of lithic and 
al 1972) are comparable to the ages ' iqn 

" , . ,. afo fbat t he flux of foramimfera > 150 

foram fluxes in DSDP site 609 also indicate that 

foram nuxes foram zones 

p dropped markedly during formabon of the first 

(Table 2). Further evidence of changes in surface water comes 

w r^m ” i — « 

sharp decreases in 51*0 of -1 per mil (Fig. 2). As th 

^ that sea surface salinities dropped 

dictate cold conditions, these results suggest th 

during deposition of all three layers. 

Based on die evidence summarized above, there can be no question 

that the carbonate-rich IRD in the four Heinrich layers has a very unusual 

origin. Any explanation based solely on North Atlantic oceanic arcu a 
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would seem difficult to reconcile with the presence of the carbonate-rich 
layers in the Labrador Sea. The evidence in hand leads us to propose a 
tentative scenario for their origin that incorporates both oceanic circulation 
and events in continental ice sheets. During each of the periods of cool 
surface water and reduced productivity coinciding with HI, 2, 4 and 5, at least 
parts of the ice sheets in eastern and northern Canada and possibly in 
northwestern Greenland as well advanced, increasing the length of the ice 
fronts in contact with the ocean. This increased the numbers of icebergs 
breaking away from ice streams that had flowed over carbonate bedrock, 
leading to the release of large amounts of carbonate-bearing glacier ice into 
the Labrador Sea. Because of the cold sea surface temperatures, melting of 
this ice was slowed and the icebergs drifted into the North Atlantic. Rapid 
transit of large amounts of this ice eastward led to rapid deposition of 
carbonate-rich IRD in a Wide belt reaching completely across the ocean. We 
have no ready explanation for the lack of abundant detrital carbonate in H3 
and H6. It could be due to less extensive advances of the Canadian-northwest 
Greenland ice fronts or to somewhat warmer sea surface temperatures 
melting away the sediment-rich basal layers of icebergs near their sources. 

Massive discharges of glacial ice might also explain the sudden 
lowering of surface salinities implied by the Sl»0 data. Hie huge amount of 
ice drifting across the North Atlantic must have produced a large volume of 
meltwater. Assuming a 8«0 of -35 X. for gladal ice (Dansgarrd et al. 1982), 
the - 1 %. S'SO change for the first three layers requires mixing only about 1 
part iceberg meltwater with about 30 parts seawater. If this is correct, large 
amounts of glacial ice must have been drifting across the North Atlantic at 
the time of H3. It is interesting that the salinity drop corresponding to the - 1 
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%, 5180 decrease in all three layers is probably enough to shut down the 

North Atlantic's thermohaline circulation. 

An important implication of the carbonate KD layers is that they 

the marine imprint of unidentified or poorly documented advances of ice 
sheets in eastern North America. What might have caused the ,ce fron 
advance is an intriguing question. If atmospheric cooling accompanied the 
shift of polar water southward, the ice advances may have been 
falling temperatures. On the other hand, as suggested by Broecker e . . 

th e lack of evidence of marked atmospheric cooling in Greenland ice cores 
the time scaie of the Heinrich layers leads to a second possibility stochastic 
surging of parts of the ice sheets. In either case, short episodes of cooling o 
N 1 Atlantic surface water appear to have been a key condition for the 
massive transit of icebergs far into the eastern North Atlantic and rapid 
deposition of carbonate-rich IRD. 
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TASK A: CLIMATE AND ATMOSPHERIC MODELING STUDIES 

riimata Model Develo p ro jeMandJ^mUcatiOM 

The research conducted during the past year in the climate and atmospheric modeling 
programs has concentrated on the development of appropriate atmospheric and upper ocean models, 
and preliminary applications of these models. Principal models are a one-dimensional radiative- 
convective model, a three-dimensional global climate model, and an upper ocean model. Principal 
applications have been the study of the impact of C0 2 , aerosols and the solar ’constant’ on climate. 

Progress has been made in the 3-D model development towards physically realistic treatment 
of these processes. In particular, a boundary layer, land surface and convection schemes have now 
been incorporated. The model containing these improvements will be run and evaluated. 

A version of the climate model has been created which follows the isotopes of water and 
sources of water (or colored water) throughout the planet. Each isotope or colored water source is 
a fraction of the climate model’s water. It participates in condensation and surface evaporation at 
different fractionation rates and is transported by the dynamics. A major benefit of this project has 
been to improve the programming techniques and physical simulation of the water vapor budget of 
the climate model. Applications include simulations of deuterium and oxygen- 18 for both current 
climate and 18,000 years ago, the source of precipitation in each grid box in the North Hemisphere, 
and a stratospheric tritium experiment to simulate the atomic testing of the 1950s and 60s (Koster el 

aL, 1990). 

During the past year, papers have been published, which investigate the impact of altered 
ocean heat transport on climate (Rind and Chandler, 1991), and the likelihood of future drought 

caused by the projected increase in temperature (Rind et al., 1990). 

Modeling of the climate and vegetation change of the last 30,000 years, and of the Little Ice 
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Age , has begun with the assistance of Rick Healy and initial attempts to compile appropriate 
boundary conditions for GCM runs. 

Studies of the sensitivity of North Pacific sea surface temperatures (SST) to global climate 
change continue. A paper in preparation will detail the results of the North Pacific cooling 
experiment, which has significant implications for northern hemisphere temperature depression, as 

well as ice sheet growth. 

An experiment for natural ocean variability on the decadal to millennium time scale has been 
conducted with a low-resolution ocean general circulation model. The model comprises a sector of 
120 degrees longitude extent, with a geometry similar to that of the Atlantic ocean. The model is 
initially run to equilibrium for 2000 years, forced with the zonal average of climatological wind stress, 
surface temperature, and salinity fluxes, parameterized in terms of a relaxation term to the zonally- 
averaged observed climatology. During the second part of the run, the surface salinity flux is 
imposed as the value obtained from the last 500 years of the first part of the run. During 2000 years 
with salinity flux forcing, the model shows considerable variability in the decadal to millennium time 
scale. The evolution of the meridional mass transport, salinity, temperature, heat and salt transport 
have been analyzed. The empirical orthogonal functions for these variables have been obtained for 
> century time-scales. A transition to one-cell pole-to-pole circulation state is observed, with 
circulation anomalies in the deep ocean appearing in the first empirical orthogonal function. The 
higher modes show variability in upper levels. Dr. Z. Garraffo (Associate Research Scientist, 
Columbia University), has done this work in collaboration with Dr. Inez Fung (GISS), using the 
Goddard Space Flight Center Cray computer. 

Work begun this past year with Dr. Silvia Garzoli (Lamont- Doherty Geological Observatory), 
on the analysis of Semter and Chervin’s global ocean general circulation model has been completed 
and published. The model output has been completed and published. The model output has been 
compared with data from regions of western boundary currents, in order to establish similarities and 
differences with observations, that will be taken into account for further use of the model results. 
The analysis contains a comparison of upper model temperature fields with satellite-derived sea 
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surface temperatures. Empirical orthogonal functions are obtained for both model and observattons. 
in order to reduce the dimensionality of the system; similarities and differences between model and 
observations are then discussed. In addition, statistical distribution functions of the observed thermal 
front me compared with model frontal positions. A comparison between upper ocean eddy kinetic 
energy predicted b, the model and fieids derived from GEOSTAT altimeter data and from drifting 
buoys has been made. The analysis was done fo, two western boundary current regions; the western 
South Atlantic and the Kuroshio-Oyashio system. In the analysed regions, maximum temperature 
differences between model and observed frontal positions were found to be in the range 2 to 4 
degrees. The eddy life cycle is slower in die model than in the observations. The model variability 
is weaker than the observed one; this difference is more pronounced in the South West Atlantic Ocean. 
Two papers from this study are in press in the Journal of Geophysical Research. 


paleocUmate -SmdigS 

Analysis of late-glacial pollen, macrofossil, acceleration mass spectrometry (AMS) controls 
and fish-scale stratigraphy from ,h. northeastern U S. has been concluded. The three cores analyzed 
indicate a sudden warming a, 12,500 ,r BP. followed by a Younger Dryas cooling between 1 1 ,000 and 
10,000 yr BP. Subsenuen, warming was also rapid, occurring in less than 10O years, and a. least three 
boreal species disappeared locally. A paper based on these findings is in preparation. Pollen and 

macrofossil research on two Vermont cores is still underway. 

Progress on research from the southeastern United States continues. Margaret Kneller (GRA. 
Columbia University) has completed the pollen and macrofossil stratigraphy from a core from 
Brown's Pond. Virginia, which contains a 17,000 yr record. She has also completed a cross-basin 
transect, comparing stratigraphy and loss-on-ignition in 5 cores. She returned to the field this past 
summer and retrieved another core in order to examine the late-glacial interval a. higher resolution. 

Research on the 18,000 year record from a Georgia core has also begun. 

Studies of climate change at high latitudes has continued, with field work in Kod.ak Island 
providing a detailed examination of the glacial to interglacial changes in stratigraphy, tephras, pollen. 
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macrofossils. an, glacial history. A P— V PoUen .a... ft- Kodiah Islam, 
climatic oscillations during late-glacial time. A. part of the invest.gat.on, .t »>« • 
whether or not this reversal is correlative with the Younger Dryas cool.ng. 

Field research in Russia (July, 1991) resulted in the acquisition of four peat cores rom . 
which me being examined for vegetational and climatic change. AMS dates on »■ > “ 

approximately 9000 yr BP, near the beginning of the Holocene Epoch. Collaboration w.th Dr. An 
Andreev (Institute of Ceography, Moscow, will result in production of a global map shown, 

initiation of high latitude peatland. 

Ongoing studies of the origin of pea, land formation in Alasha have provided AMS dates on 
macrofossils from permafrost pea, cores. The pea, formation ranges in age from 16,000 ,0 g.OOO yr 

BP. 


Climate C hanES_ApJPUoatiailS. . . !f 

A Coastal Hazards Data Base for the Southeast U.S. has been created to categorize and identify 

shorelines a, high risk to future sea level rise, caused by climate warming. The data base consists of 

clusters, which, in order of importance, represent permanent inundation. ep.sod.c 
erosion potential, respectively. Each data variable was assigned to one of 3 rish Casses, and these 
were combined into a Coastal Vulnerability Index (CV.) that provides a relative measure of 
faced b, any given segment of the shoreline. Among the highest rish shorelines in «h K reg,on me 
Cane Hatteras, North Carolina and the coast of Louisiana. A, Cap. Ha, .eras, the episodic .nun anon 
variables (most,, climate variables) contribute a significant,, larger proportion ,0 the torn 
whereas a, Caminada Pass, the permanent inundation factors (due to the anomalously h.gh ) 
assumes a greater share of the total score. Sea level rise (SLR) scenario,, modified from , he 1 
and 1992 findings, for local subsidence effects, were applied to six , est-si.es in the US w . 
have been selected on the basis of the CV, rating. These SLR scenarios were utilized to calcu^te the 
extent of inundation under various shoreline response models, as documented more fully .n 
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et al. (1992). 

A project was initiated to assess the importance of methane hydrate in permafrost and deep- 
sea sediments as a future potential source of atmospheric methane, under warmer climate conditions. 
In the initial phase of this study, estimates have been made of methane hydrate reservoirs m the 
world’s oceans, based on two proposed models of hydrate formation: a) in-situ bacterial formation 
and b) pore fluid-migration. While both models predict similar methane volumes, their distribution 
of gas hydrate with respect to the sub-seafloor surface differs significantly in terms of the thermal 
response to climate warming. Since the thermal pulse from increased ocean temperatures will travel 
downward into the seafloor sediment, hydrates distributed uniformly with depth from the ocean 
bottom, as in the first case, will potentially release more methane within the next 100 years, than 
hydrates concentrated at depth, as predicted by the second model. Further work is in progress to 
assess the relative merits of these two models and to calculate the thermal penetration at depth, for 
plausible scenarios of global temperature warming within the next 100 years. 

SAGE II 

In order to explore the solar cycle/tropospheric connection, several preliminary experiments 
have been conducted to assess the likelihood of producing a realistic response. Considering that the 
effect of solar activity variations is apparently observed more clearly when modulated by the phase 
of the Quasi-Biennial Oscillation (QBO), in the first experiment, a QBO was induced in the model, 
by using an additional momentum forcing in the lower stratosphere to produce either strong tropical 
east winds or west winds (the model under normal circumstances does not generate a QBO, perhaps 
because of insufficient resolution). Then a set of initial conditions (outside the tropics) which had 
produced a stratospheric warming in the control run was reintroduced. The model generated a strong 
warming in the QBO east phase, but not in the west phase. This was the result of increased upward 
wave activity flux from mid-latitudes during the QBO-east. The tendency for warming mainly 
during the east phase is in agreement with observations. The experiments suggest the model is able 
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to differentiate dynamic events 


in its different QBO phases despite the artificiality of the QBO 


generation 


the hex. experiment, using the control run. the UV radiation below of 0.3u was 1 

by 50%. This is greatly in excess of observed values (which are still somewhat uncertain) but tn line 
with a general modeling approach: firs, use strong forcing to elucidate mechanisms, and see if there 
is an, response, then reduce the forcing bach to realistic magnitude,. With the increased UV heating 
of the ozone layer in the tropical regions receiving solar insolation, the west wind, in the Northern 
Hemisphere increased (consistent with the therma. wind relationship,. The stronger zona, winds 
reduced the vertica, propagation of wave activity and reduced the amplitudes of planetary waves. 

TASK B: CLIMATE APPLICATIONS OF EARTH AND PLANETARY OBSERVATIONS 
Cloud Clim atology: 

During the pas. year, the study of tropica, convection using 1SCCP data has continued. 
Examination of changes in convection both during El Nino events and between non-El Nino years 
shows .ha, the nature of the shift in convective properties associated with longer time scales ts no. 
the same as the rapid variations that occur on daily to monthly time scales. Moreover, the E. Nino- 
related changes appear similar to interannual changes. An analysis of possible feedback on tropica, 
sea surface temperatures (SST) was performed to examine a hypothesis that tropical SST ts controlled 
b, convection through a local cloud-radiative feedback. The msults show that the ne, effect of 
changing convective clouds on surface radiative heating is small compared with changes tn 
evaporative cooling and tha, neither change is large enough to completely control the SST changes. 
This confirms other studies that show tha, El Nino SST changes me con, roiled more by ocean 
dynamics than by changes in the surface energy budge,. These results have been published. 

A global survey of cloud particles sizes for low-level. linuid water clouds was conducted for 
the four seasons, and for morning and afternoon, using multi-spectral radiance dam from two polar 
orbiting weather satellite,. Although validation dam are scarce, agreement with in-si, u measurements 
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from FIRE (First ISCCP Regional Experiment) was found to be excellent, for two cases. The results 
of this survey confirmed eariier aircraft results in showing tha, cloud droplet sixes are systematica.., 
larger in marine clouds as compared to land clouds. Moreover, a systematic hemispheric difference 
was found in the droplet sixes, which appears to he consistent with that predicted by the different 
abundances of tropospheric aerosols. A significant diurnal variation of cloud droplet size, increasing 
from morning to afternoon, is also found to occur in marine and land clouds at low latitudes. Thts 
study shows tha. the relationship of cloud particle size, optical thickness and water content is more 
variable than commonly assumed in most climate models. A paper btrned on this Ph.D. thesis research 

of a recently-graduated GRA is in preparation. 

A study of the variations of the optical thickness of low level clouds a. low and middle 

.atitudes has been extended to a combination of satellite and surface weather observations. The 
earlier results, based on satellite data, have been submitted for publication. They showed that there 
is change in the behavior of low-level cloud optica, thickness with temperature. The optical 
thicknesses of "cold- clouds (T < -5'C) increase with increasing temperatures, a rate that is 
proportional to the increase of the adiabatic cloud water amount in clouds. On the other hand, the 
optica, thicknesses of -warm- clouds decrease with increasing temperatures. If this behavior aiso 
occurs during a climate warming event, then the global cloud radiative feedback is positive. The 
cloud-radiative feedback also varies with latitude in such a fashion as to nearly eliminate the 
amplification of climate change with latitude that is found in most climate models, but m no. 
observed. The extended study addressed the question of whether these observed changes in cloud 
op, ica. thickness were due tochanges in cloud water con, entorwhether Urey were caused by changes 
in cloud droplet size or vertical extent. Satellite measurements of cloud particle sizes indicate tha, 
the changes are too small to explain the observed changes in optical thickness. Combining rawmsonde 
and satellite data shows that vertical extent changes also do not explain the observed optic 
changes. Moreover, the hypothesis that the change in behavior of cloud water content is as 
with a change in precipitation efficiency with temperature is supported by the changing frequency 
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of precipitation occurrence with temperature found from surface observations. These results 
represent the Ph D. thesis research of a second GRA that was successfully completed this year, and 

form the basis of another paper, in preparation. 

New convection and cloud parameterizations for the GISS climate GCM continue to be tested 

by comparisons against the ISCCP data and by performing climate sensitivity studies with the climate 
model. One notable success is the clear improvement of the model cloud changes associated with an 
El Nino even, (where the model is forced by the observed changes in SST) when die new convection 

and cloud parameterizations are used. 

Comparisons of cloud water contents inferred from ISCCP cloud properties and from 
microwave measurements show excellent agreement for liquid water clouds (the microwave 
measurements are no. sensitive to the ice phase), except for very low values. The microwave 
measurements have difficulty detecting low water contend because of uncertainties in water vapor 
effects; hence the optical measurements of ISCCP are more sensitive mid 'see- about twice as much 
cloud as the microwave measurements. The high resolution measurements of ISCCP can also be used 
to evaluate other sources of error for the microwave measurements. These resold are being prepared 

for publication. This work represend the thesis research of a third GRA. 

A climatology of cloud layer thicknesses has been obtained that is based on 14 years of 
rawinsonde measuremend of temperature and humidity profiles from 60 northern hemisphere sites. 
This represend a unique dataset. Comparison of cloud statistics in this dams., with the full cloud 
climatology from surface observations indicates that the cloud layer climatology is generally weU 
sampled; however, it is composed almost exclusively of cases with a single cloud layer. Analysis of 
the resold indicates little diurnal or seasonal variation in cloud layer thicknesses but a systematic 
increase of thickness with latitude and altitude. However, a more thorough analysis is needed, 
because the resold clearly depend on cloud types; the actual layer thickness distributions within each 
altitude range are multi-modal. The preliminary analysis of these resold is being prepared for 
publication, based on the thesis research of a fourth GRA. 
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Three significant improvements have been made to the procedure for calculating global, total 
solar and thermal infrared fluxes at the top of the atmosphere and at the surface from the 1SCCP 
dataset. The most important change is in the prescription of the surface albedo from the ISCCP 
visible reflectances: the estimated near-IR albedos for land surfaces proved to be too large when 
compared with ERBE observations have been calculated. A new climatology of tropospheric aerosol, 
with more realistic geographic variations, has been added to the radiative model. A more accurate 
parameterization of the angular dependence of radiative transfer through clouds has been 
implemented in the radiative model. All of these changes have improved the validation comparisons. 
Validation studies and sensitivity tests have been completed and the results are being prepared for 
publication. Routine calculations, to cover several years of data, have commenced. These results will 
be used to study cloud radiative feedbacks on synoptic, seasonal and longer-term variations. 

A simpler method for calculating surface solar irradiance was developed and described in a 
paper (Bishop and Rossow 1991). A 7.5 year climatology of daily surface solar irradiances has been 
produced and is being analyzed for ocean biospheric studies. 


Tracking of polarization features above the visible cloud level on Venus yields a mean zonal 
velocity of -23.5 m/s, which is much slower than derived from tracking of UV features near the 
cloud top. Preliminary estimates of haze optical depth reveal two distinct regimes: one of high haze 
optical depth and high variability in polar regions, and the other of lower optical depth and variability 

at lower latitudes. 
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The Goddard Institute for Space Studies (GISS) provides funding for Lamont-Doherty 
Geological Observatory (LDGO) to support graduate student research in geochemistry with close ties 
to subjects under study by Dr. James Hansen and the Climate Modeling group at GISS. In most cases, 
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these projects involve use of the GISS GCM. This progress report briefly outlines the activities of 
eight current graduate students who are working with Professor Wallace S. Broecker, Geochemistry. 

at LDGO on projects related to global change. 

Fritz Zaucker (GRA, LDGO) is expecting to defend his thesis on the water vapor transport 

through the atmosphere and to impact on ocean circulation, in dr. Fall. 1992. He Intends to accept 
a NASA post-doctoral fellowship at GISS. 

Kevin Harrison (GRA, LDGO) is starting his final year at Columbia University. His thesis 
is directed toward estimating the eaten, ,0 which CO, -induced growth enhancement has increased 
the inventories of carbon in soils, litter, and vegetation. A paper describing his work to date is 

included in Appendix A. 

Jeff Severinghaus (GRA, LDGO) is just beginning his Ph.D. research. He is doing trace gas 
and isotope measurements in Biosphere 2 in order to document the utility of research relating 
materials from closed ecosystems to our understanding of terrestrial biogeochemistry, 
project is ,0 determine the ratio of ACO, to AO, and A S ”C ,0 ACO, throughout the diurnal cycle. 
The pCO, in Biosphere 2 current goes from 1300 ppm at day’s beginning to 500 ppm a, day’s end. 
His second project will be to determine the Dole effect (ozygen isotope fractionation during 

repiration). 

Roberto Gwiazda (GRA, LDGO) is in the second year of Ph.D. research. He rs trying to 
assess the role of soil CO, in the rate of silicate weathering. To date, he has been doing pilot studies 
on small basins in the vicinity of Lamont (soil CO,, ground water CO, and stream alkalinity). 

Lucy Charles (GRA. LDGO) is a first year student. She has been working with Drs. Dave 
Archer and Taro Takahashi on mixed layer ocean chemistry modeling and with Prof. W.S. Broecker 
on evaluating the estimates of the northern to southern hemisphere transport of carbon in the forms 

of CH 4 and CO. 

Jerry McManus (GRA. LDGO) is a second year student. He is working on the ice rafting 
record in northern Atlantic sediments. In particular, he is trying ,0 determine the causes and effects 
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of large armadas of icebergs released into the Atlantic by surges of the eastern margin of the 
Laurentian ice sheet (see Appendix B). 

Abhijit Sanyal (GRA, LDGO) is a first year Ph.D. student. He is attempting to document the 
extent of chemical erosion in the deep tropical Pacific during Holocene time. This study has 
important implications to the debate about the role of the CaCOj cycle in glacial to interglacial 
atmospheric CO a content changes. 

Rachel Oxburgh (GRA, LDGO) is a second year Ph.D. student. She is trying to reconstruct 
oceanic osmium isotope content changes over the last 60 million years through measurements on deep- 
sea sediments. The purpose is to determine how global chemical weathering rates changed forty 
million years ago when the Himalayas were created through the collision between the Indian and Asia 
land masses. 


(1991-1992 and in press) 
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